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ABSTRACT 

A survey of Type II supernovae explosion models has been carried out to determine how their 
light curves and spectra vary with their mass, metallicity and explosion energy. The presupernova 
models are taken from a recent survey of massive stellar evolution at solar metallicity supplemented 
by new calculations at sub-solar metallicity. Explosions are simulated by the motion of a piston near 
the edge of the iron core and the resulting light curves and spectra are calculated using full multi- 
wavelength radiation transport. Formulae are developed that describe approximately how the model 
observables (light curve luminosity and duration) scale with the progenitor mass, explosion energy, 
and radioactive nucleosynthesis. Comparison with observational data shows that the explosion energy 
of typical supernovae (as measured by kinetic energy at infinity) varies by nearly an order of magnitude 
- from 0.5 to 4.0 x 10 51 ergs, with a typical value of ~ 0.9 x 10 51 ergs. Despite the large variation, 
the models exhibit a tight relationship between luminosity and expansion velocity, similar to that 
previously employed empirically to make SNe IIP standardized candles. This relation is explained by 
the simple behavior of hydrogen recombination in the supernova envelope, but we find a sensitivity to 
progenitor metallicity and mass that could lead to systematic errors. Additional correlations between 
light curve luminosity, duration, and color might enable the use of SNe IIP to obtain distances accurate 
to ^20% using only photometric data. 

Subject headings: distance scale - radiative transfer - supernovae: general 



1. INTRODUCTION 

Type II supernovae (SNe II) result from the explosion 
of massive stars that have retained their hydrogen enve- 
lope until their cores collapse to neutron stars or black 
holes. The most common events, the Type II plateau 
supernovae, have a distinctive light curve, maintaining 
a nearly constant luminosity for ~ 100 days, then sud- 
denly dropping off. Upcoming synoptic surveys should 
discover millions of these events out to redshifts of a few. 
These observations will probe massive stellar evolution 
in a broad range of galactic environments, and may also 
be used to measure cosmological distances and test host 
galaxy dust properties. 

SNe II are diverse transients in terms of their luminosi- 
ties, durations, and expansion speeds. By modeling the 
observed light curves and spectra, one can constrain the 
physical properties of the explosion such as the ejected 
mass, explosion energy, and pr esupernova radius. In a pi- 
oneering study, iHamuvl (|2003f ) analyzed 16 observed su- 
pernovae and derived masses between 15 and 50 M Q and 
presupernova radii ranging from 70 to 600 R© . Unfortu- 
nately these inferred values seem implausible, especially 
the high masses. Direct observations of the progenitors of 
nearby SNe II fro m pre-explosion ima ges indicate masses 
of only 8—15 M Q (jSmartt et al.ll2008h . while modern stel- 
lar ev olution models predi c t a mass range of ~ 12— 25 M 
(e.g.. IWooslev et alJl2002t IWooslev fe Hege3l2007h . The 
overly large values inferred by IHamuvl (|2003l ) likely re- 
flect deficiencies in the theoretical models used in the 
light curve analysis (jLitvinova fe Nadczhinl ll985f h More 
recent models tailored to individual e vents have returned 
more reasonable numbers (e.g. lUtrobin fe Ch ugai 2008; 
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lBaklanovet~aIll2005l) . 

Observations of SNe II are also used to mea- 
sure cosmological distances. In terms of bright- 
ness alone, they are poor standard candles with lu- 
minosities varying by more than an order of mag- 
nitude, but various methods can be used to stan- 
dardized them. Most previous studies were of the 
B aade-W essellink type, e.g. , the expanding photo s phere 
method dKirshner fc Kwanl 119741 lEastman et all 119961: 
i Dessart fc Hillierll2005at iDessart et al.ll2008t Pones et al.1 
2009) and the spectral-fitting expanding atmospher e 
method fSEAM. lMitchell et alj|20"0"a iBaron et al.ll2004h . 
Both approac hes require detailed atm ospheric modeling. 
More recently, Hamuy fc Pintol (|2002l ) suggested a much 
simpler approach using an empirical correlation between 
plateau luminosity and expansion velocity (as measured 
from the Doppler shift of spectral lines). This standard- 
ized candle method has been exte nded out to moder- 
ate redshifts with pro mising results (jNugent et al.l 120061 : 
iPoznanski et al.ll2008f ). However, the physical underpin- 
nings of the method and its potential vulnerability to 
systematic error have not yet been fully explored. 

Here we present a new survey of SNe IIP models based 
on the one-dimensional explosion of realistic progenitor 
star models of various masses, metallicities, and explo- 
sion energies. We calculate the broadband light curves 
and the detailed spectral evolution using a code that in- 
cludes a full multi-wavelength solution to the radiation 
transport problem. The models allow us to explore the 
light curves dependence on model parameters, which we 
compare to analytic expectations. We also reproduce and 
illuminate the physics leading to the standardized candle 
relation. 
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2. ANALYTICAL SCALINGS 

Some insight into the light curves of SNe IIP can be 
gained from analytical scalings which express how the ob- 
servables - luminosity, L sn , light curve duration, t sn , and 
expansion velocity, v sn ~ depend on the basic supernova 
parameters - the explosion energy, E, ejecta mass, M e j 
and presupernova radius, Ro- Although the analytic light 
curves themselves are only a pproximate, the scaling re- 
lations can b e quite accurate (|ArnettHl98Ct lChugailll991l : 
lPopov|[l993ft . As a guide to the model results below, 
we rederive here some basic results using simple physical 
arguments. 

Although the mechanism of core-collapse supernova ex- 
plosions is complicated, the end result is the deposition 
of order E ~ 10 51 ergs = 1 B near the center of the star, 
and subsequently, the propagation of a blast wave that 
heats and ejects the stellar envelope. At the time the 
shock reaches the surface (hours to days), the explosion 
energy is roughly equally divided between internal and 
kinetic energy. In the hydrogen envelope, radiation en- 
ergy dominates the thermal energy of ions and electrons 
by several orders of magnitude. 

In the subsequent expansion, the internal energy is 
mostly converted into kinetic. The ejecta is optically 
thick to electron scattering, and the adiabatic condition 
gives E int (t) = E a (R(t) / Ro)' 1 , where E w E/2 is the 
initial internal energy, and R, the radius. After many 
doubling times, the ejecta reaches a phase of free homol- 
ogous expansion, where the velocity of a fluid element is 
proportional to radius R = vt. The final velocity of the 
ejecta is of order 



(2E/M cj ) 1/2 = 3 x 10 8 M W E 51 cm s~ 



(1) 



where M 10 = M cj /10 M Q , E 51 = E/10 51 ergs. For 
homologous expansion, the internal energy evolves as 
Ei n t{t) = Eo(t/t e )~ l where t e = Rq/v sii is the expansion 
time. Dimensionally, the luminosity of the light curve 
will then be 



Lay, 



Eint(t S n) _ Eot e 



t 2 



(2) 



where t sn is the appropriate timescale for the duration 
of the light curve. We determine i sn and t sn below for 
three different scenarios. 

First, if the opacity, k, in the supernova envelope were 
a constant, the time scale of the light curve would be 
set by the effective diffusion time. Given the mean free 
path A p = (up)' 1 and the optical depth of the ejecta 
r = R/Xp, the diffusion time is 



is 



2 A p 



R 2 k P 



(3) 



Over time, the supernova radius increases and the den- 
sity decreases due to the outward expansion. Using 
the characteristic values R(t sn ) — v sn t sn and p(t sn ) ~ 
M C j/R(t sn ) 3 in Eq. 02 one can solve for i sn . 



t sn OC E 



L sn cx E M^R k-\ 



(4) 



where the scaling for L sn was d etermined by u sing t sn in 
Eq. [2] These are the scalings of lArnetti 0980) 



The results Eq. 31 though commonly applied, are not 
quite adequate for SNe IIP because the assumption of 
constan t opacity neglects the important effects of ion- 
ization (jGrassberg et al.lll971[ K Once the outer layers of 
ejecta cool below T\ w 6000 K, hydrogen recombines and 
the electron scattering opacity drops by several orders of 
magnitude. A sharp ionization front develops - ionized 
material inside the front is opaque, while neutral ma- 
terial above the front is transparent. The photon mean 
free path in the ionized matter (A m f p ~ 10 10 cm) is much 
smaller than the ejecta radius (R = v s t « 10 15 cm), so 
the supernova photosphere can be considered nearly co- 
incident with the ionization front. 

As radiation escapes and cools the photosphere, the 
ionization front recedes inward in Lagrangian coordi- 
nates, in what is called a recombination wave. The pro- 
gressive elimination of electron scattering opacity allows 
for a more rapid release of the internal energy. Once the 
ionization front reaches the base of the hydrogen enve- 
lope, and the internal energy has been largely depleted, 
the light curve drops off sharply, ending the "plateau" . 
Any subsequent luminosity must be powered by the de- 
cay of radioactive elements synthesized in the explosion 
(the light curve "tail"). 

To account for hydrogen recombination while ignoring , 
for the mome nt, radiative diffusion (as in lWoosievlll988l ; 
IChugailH99l . we use the fact that the photosphere is 
fixed near the ionization temperature T\ and radiates a 
luminosity L p given by 



L p = 47ri? 2 crT 4 



(5) 



The light curve timescale t sn is then given by how long 
it takes the photosphere to radiate away all of the (adi- 
abatically degraded) internal energy - in other words, 
Eo(t e /t m ). This, (along with R = v s t sn ) gives 



Lpt sn 



t sn ocE-^Ml{ & Rl^T-\ 
L sn oc E^M^R^T 2 . 



(6) 



Unlike Eqs. [H these scalings show no dependence on k. 
Effectively, the assumption is that the ejecta is infinitely 
opaque below the recombination front, and fully trans- 
parent above. We will sec in the numerical models that 
this assumption is not totally correct, and hence diffusion 
below the photosphere is important. 

Finally, to derive scalings which include both the ef- 
fects of rad iative diffusion and recombination (as in 
lPopoviri993h . we return to the diffusion time equation 
of Eq. [3l but now realize that the radius of the opaque 
debris changes over time, not only due to the outward 
expansion but also from the inward propagation of the 
recombination front. This photospheric radius is deter- 
mined from Eq. [5] 



R 2 = Lp 



Ent. 



' e 



4ttoT? t 2 sn 4TT,aT? 



(7) 



where in the last equality we used Eq. [2] to rewrite 
L p = L sn . Plugging this expression for the radius into 
the numerator of Eq. [3] gives 

t sn oc S-^M, 1 / 2 ^ V/ 6 T- 2/3 

^5/6 A , f -l/2 p 2/3 -l/3 T 4/3 



L sn (X ^"M^'^'^^X 



(8) 
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TABLE 1 

Progenitor Star Models 





Z 


Mi 


Ro 


M Fc 




12 


1.0 


10.9 


625 


1.365 


0.30 


15 


1.0 


12.8 


812 


1.482 


0.33 


15 


0.1 


13.3 


632 


1.462 


0.33 


20 


1.0 


15.9 


1044 


1.540 


0.38 


25 


1.0 


15.8 


1349 


1.590 


0.45 



These scalings are identical to those found by iPopovl 
(1993) in a more involved analysis. 

So far, the light curves described have only accounted 
for shock deposited energy. Radioactive 56 Ni synthesized 
in the explosion introduces an additional energy source 
concentrated near the center of the debris. The heating 
from radioactive decay helps maintain the ionization of 
the debris and so extends the duration of the plateau. 
To incorporate this effect into the scalings, we generalize 
the expression for the internal energy 



Eintit) — E — 



(9) 



where E ni w 0.6 x 10 50 Af Ni ergs, E co w 1.2 x 10 50 Af Ni ergs 
are the total energy released from 56 Ni and 56 Co decay 
(with M-Ni in units of M Q ) and i n ; w 8.8, i co « 113 days 
are the lifetimes. This correction for radioactivity can be 
written E int (t) = E (t e /t sn )f lad with 



fr 



ad 



1 + 0.26 



E5I te 



(10) 



where E51 = i?/10 51 ergs. Following the same arguments 
leading to Eq. [3J we see that the plateau timescale scales 

as t sn oc /j-Yd- F° r example, a 56 Ni mass of 0.1 Mq should 
extend the plateau by ~ 24% for E51 = 1. Although one 
might anticipate a change to L sn as well, the models show 
that the decay energy does not typically have enough 
time to diffuse out and affect the plateau luminosity (see 
Section O . 

Below we compare the analytic scalings to our numer- 
ical simulations and show that the simple relations, par- 
ticularly Eqs. [5J agree quite well. The models allow us 
to refine the exponents and calibrate the numerical con- 
stants in front. 

3. PROGENITOR MODELS 

For our numerical models, we consider stellar progen- 
itors with main sequence masses in the range Mj = 
12 — 25 M Q , the range expected to produce most of the 
observed events. Properties of the presupernova models 
are summarized in Table 1 which gives the zero age main 
sequence and presupernova masses in solar masses (Mi 
and M{), the presupernova radius in solar radii (Ro), the 
iron core mass in solar masses (Aipc), and the surface 
helium mass fraction. All models were computed with 
the Kepler code, which follows stellar evolution includ- 
ing the most up-to-date opacities, prescriptions for mass 
loss, and nuclear reaction rates dRauscher et alJ 120021 : 
IWooslev et aL|[2002l: iWooslev fc Hegerll2007l ). Stars with 
larger initial masses experience more mass loss, especially 
during the red giant phase, and this narrows the range of 
final masses to M f = 10.9 - 15.8 M Q . For M\ > 20 M Q , 




Fig. 1. — Abundance distribution of model M15-E.12 after explo- 
sion assuming either no mixing (top panel) or the mixing applied 
in this paper (bottom panel). The inner hashed region form the 
remnant.. 

the presupernova mass declines with increasing M\. More 
massive stars do, however, maintain significantly larger 
radii at the time of explosion. 

The helium mass fraction in the stellar envelope is 
also a function of M\ (Table 1) varying from 30% for 
Mi = 12 M to 45% for Af ; = 25 M . This variation 
is due to mass loss and convective dredge up from the 
helium core, which are greater in more massive stars. As 
the envelope helium abundance affects both the electron 
scattering opacity and the recombination temperature of 
the ejecta, we will find it has a significant effect on the 
light curves of SNe IIP. 

While most stars in our survey have solar abundances, 
a lower metallicity (Z = 0.1 solar) 15 M Q model was also 
included. The chief effect of the lower metallicity was a 
smaller presupernova radius and less total mass loss. 



4.1. 



4. EXPLOSIONS AND NUCLEOSYNTHESIS 
Mass Cut, Fallback, and the Production of 56 Ni 



The explosion of each model was simulated by moving 
a piston outwa rd from an inner boundary at mass co- 
ordinate Af pist (fWooslev fc Weaverl ["19951 ; IWooslev et al.1 



2002f) . typically taken to be the outer edge of the iron 
core, and following the subsequent hydrodynamics as- 
suming radial symmetry. Each star was exploded several 
times to obtain variable kinetic energies at infinity within 
the set of approximately 0.3, 0.6, 1.2, 2.4, and 4.8 x 10 51 
erg. The results are summarized in Table 2. Polarization 
observations of SNe IIP suggest that the hydrogen en- 
velopes are indeed spherically symmetric, although the 
cores may appear aspherical , perhaps due to an asym- 
metric explosion mechanism (jLeonard et a l. 2006, 2001). 
Any asymmetry in the shock wave, however, is likely 
smoothed out by propagating through the large hydrogen 
envelope. 

Explosive nucleosynthesis was calculated us i ng th e 
same code and physics as in IWooslev fc Heeen ([2007h . 
While the mass of 56 Ni that is synthesized is numerically 
well determined by this procedure, it is an overestimate 
for two reasons. First, situating the piston at the edge of 
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the iron core, the deepest it can possibly be without vi- 
olating nucleosynthetic constraints on the iron isotopes, 
overestimates both the density and mass close to the ex- 
plosion. It is difficult to launch a successful explosion in 
the face of such high accretion and a more reasonable lo- 
cation for the piston might be farther out near the base 
of the oxygen shell. There is a sudden increase in the en- 
tropy per baryon, S/JV^fc, to a value around 4.0 signals 
a rapid fall off in density in the presupernova star. To 
illustrate the sensitivity to the piston location, a second 
version of the 0.6 B explosion of the solar metallicity 15 
Mq star was calculated with the piston at the location 
where S/Na^ = 4.0. The 56 Ni production declined from 
0.24 M to 0.084 M Q . 

Second, the ejection of 56 Ni is sensitive to the treat- 
ment of mixing and fallback. It was assumed here that 
whatever material had positive speed at 10 6 s, the time 
of the link from the Kepler hydrodynamics code to the 
spectral synthesis code would be ejected. Since some of 
the slow moving material will fall back into the collapsed 
remnant at later times, this also overestimates the ejec- 
tion of 56 Ni. However, mixing which is largely finished 
before 10 6 s, reduces this sensitivity by taking 56 Ni that 
would have fallen in and moving it farther out in the 
ejecta (jHerant fc WooslevlH99llJoggerst et al.ll2009f ). 

Observations of the tail of the supernova light curve 
constrain the mass of 56 Ni to be in th e rang e ~0.01 - 0.1 
M Q (jArnett et all 1 19891: ISmartt et all 120081 ) with values 
closer to 0.1 Mq coming from the more massive progen- 
itors. All in all, it seems that multiplying the 56 Ni yield 
of our models by a factor of 0.25 — 0.5 is reasonable. It 
should be noted that the plateau phase of the light curve 
is insensitive to the 56 Ni production (see Section [6]). 

4.2. Mixing 

In addition to its affect on the absolute yields, hydro- 
dynamical mixing during the explosion can carry 56 Ni 
out into the hydrogen envelope and hydrogen deep into 
the core of helium and heavy elements. The early ap- 
pearance of X-rays in SN 1987A and the smoothness of 
the light curve showed that substantial mixing occurred 
- more than has been provided so far in any calculation 
of just the Rayleigh- Taylor instability. Mixing that com- 
mences with a broken sym metry in the exploding core 
itself seems to be necessary (jKifonidis et al.ll200l . Be- 
cause a large number of models needed to be studied here 
and because the degree of mixing is affected by the uncer- 
tain asymmetry of the central engine, we used a simple 
paramet ric representation of the mixi ng similar to that 
used by IPinto fc Wooslevl ([1988) and IHeger fe Wooslevl 
(2008]). A running boxcar average of width AM is moved 
through the star a total of n times until the desired mix- 
ing is obtained. The default values AM and n are 10% 
of the mass of the helium core and 4, respectively. This 
gives, for example, the mixed composition for Model 15C 
in Figure [T] 

We explored the effects of varying the degree of mixing, 
and found that it lead to only small changes at the end 
of the plateau - i.e. once the recombination wave had 
reached the inner layers of helium and heavier elements. 

5. THE CALCULATION OF LIGHT CURVES AND SPECTRA 

Several numerical studies of the light curves of SNe IIP 
have been published (e.g.. iLitvinova fe Nadezhinl 119851 : 



42.5 - 




Days since explosion 

Fig. 2. — Bolometric light curves of a standard model 
(M15_E1.2_Z1) with different amounts of 56 Ni ejected (marked on 
figure in units of Mq). The radioactive energy deposition extends 
the plateau, but has little impact on the luminosity at t < 50 days. 

I Young fc Branch! [19891 : ITltrobinl [20071 : INadvozhinl[200l . 
One common limitation of the previous studies was 
that the radiative transfer was often treated in the 
diffusion approximation or with low wavelength res- 
oluti on, although there have been a few ex ceptions 
(c.f., IBaklanov et "all [20051 : IChieffi et alJ l200l . Non- 
LTE (NLTE) radiative transfer calculations have been 
applied to the stationary spectra of SN e IIP (e.g., 
iDessart fc Hilher1l2005bl : iBaron et all [20031 ). but not, so 
far, to time-dependent light curve calculations. 

To calculate light curves and spectra of our mod- 
els, we applied a novel method which coupled a 
multi-wavelength implicit Monte Carlo radiation trans- 
port code to a 1 -dimensional hydrodyn amics solver 
ijKasen et al][200l iKasen fc Wooslevl [2009] ). The initial 
conditions of the calculation were taken from the Ke- 
pler explosion model at t — 10 6 s after explosion. At 
this time, the ejecta was largely homologous and the hy- 
drodynamics essentially unimportant. While we there- 
fore neglect the earliest part of the light curve, our main 
interest here is the plateau phase. Detailed radiative 
transfer calculations of the shock breakout phase and 
early luminosity will be discussed in a separate paper 
(| Kasen fc Wooslevll2009l) . 

In the Monte Carlo approach, the radiation field is rep- 
resented by discrete photon packets which are tracked 
through randomized scatterings and absorptions. At the 
start of the calculation, a large number (~ 10 5 ) of pack- 
ets were initiated in each zone. The energy of the packets 
was chosen so that the sum equaled the equilibrium ra- 
diation energy of the zone. The initial frequency v and 
direction vector D of each packet were sampled assuming 
that the distribution was isotropic and blackbody in the 
comoving frame. Throughout the simulation, additional 
packets were created to model gamma-rays input by the 
decay of 56 Ni and 56 Co. The transport and absorption of 
these gamma rays were likewise followed using a Monte 
Carlo approach applying the relevant opacities. 

We adopted a mixed-frame approach for the transport 
whereby the gas opacities and emissivities were calcu- 
lated in the comoving frame, while Monte Carlo packets 
were tracked in the lab frame. The relevant optical opac- 
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ities included electron scattering, bound-free, free-free, 
and bound-bound line opacity, the last treated in the ex - 
pansion opacity formalism of lEastman fc Pintol (|1993f) . 
The matter ionization and excitation state were com- 
puted assuming Saha/Boltzmann statistics at the matter 
temperature. While the code allows for non-equilibrium 
between matter and radiation temperatures, the radia- 
tion energy density is so dominant in SNe IIP that the 
two equilibrated on a short timescale. 

The scattering of photon packets was simulated by 
Lorentz transforming a packet into the comoving frame, 
preforming an isotropic scattering, and then transform- 
ing back to the lab frame. The application of the two 
Lorentz transformations changes the energy and fre- 
quency of the outgoing packet. When averaged over 
many scattering events, this effect accounts for the work 
done by the radiation field. We checked that the correct 
behavior was recovered in very optical thick regions of 
the ejecta, where the radiation energy density evolved 
with time as it should for a homologous adiabatic flow, 
e ra d oc (t/t ) -4 - 

While the properties of individual packets were sam- 
pled from continuous distributions in space, time, and 
wavelength, the grid through which they moved was dis- 
crete. In these calculations, the ejecta was divided into 
150 equally spaced radial zones. Opacities and emissiv- 
ities in each zone were further defined on a wavelength 
grid of range 1 — 25, 000 A with a constant binning of 
5 A. The physical properties of the zones (e.g., density, 
temperature, ionization state, and opacity) were updated 
on a timescale chosen much shorter than the dynamical 
timescale, with time-steps not exceeding 5x 10 4 s. Higher 
resolution tests were performed to confirm that the dis- 
cretization was adequate. 

NLTE calculations of Type II spectra show that devi- 
ations from LTE have significant effects on line profiles, 
while the continuum flux is less affected (jBaron et al.l 
U996tlDessart fc Hillierll2008f ). To estimate the potential 
effects, we computed stationary NLTE spectra on the 
plateau (day 50 after explosion) using the same code. 
A particularly relevant NLTE effect is on the Ca II IR- 
triplet, whose emission is over predicted by LTE enough 
to cause a ~ 0.1 mag increase in /-band magnitude. 
While still relatively small, this error could be signifi- 
cant when using the predicted V — I color to correct for 
dust extinction in observations. We therefore used the 
NLTE spectral results for the discussion in Section [HI 
Time-dependent NLTE effects, which are not included 
here, can stro ngly affect the emissio n in the hydrogen 
Balmer lines (|Dessart fc Hillierl 120081) and the H-alpha 
line in particular, which would modify the i?-band mag- 
nitudes. 

6. A TYPICAL MODEL 

6.1. Bolometric Light Curve 

We first focus on the properties of a typical SNe IIP 
model (M15_E1.2_Z1) which has parameters thought to 
be common: Mj = 15 M Q , E — 1.2 B, and solar metallic- 
ity. Figure [2] shows the bolometric light curve for differ- 
ent values of the ejected 56 Ni mass. Initially, the model 
luminosity decreases after shock breakout, and reaches a 
minimum around day 20. At that time, the outermost 
layers of ejecta become cool enough such that hydrogen 




Velocity (km/s) 

Fig. 3. — Evolution of the hydrogen ionization fraction over time 
for model M15_E1.2_Z1. The ejecta photosphere forms near the 
ionization front, which recedes inward in velocity (i.e., Lagrangian) 
coordinates as the photosphere radiates and cools. The plateau 
ends when the ionization front reaches the base of the hydrogen 
front. 




1 ... i ... i ... i ... i ... i ... i ... i . 
20 40 60 80 100 120 140 
days since explosion 

Fig. 4. — Photospheric velocity evolution of model M15_E1.2_Z1 
as a function of time. The solid line shows the velocity of the 
electron scattering photosphere defined at r e = 2/3. The dashed 
line shows the photospheric velocity in a model in which hydrogen 
recombination is neglected. The circle s are the obse r ved v elocities 
of the sample of SNe IIP compiled bv lNugent et al.l 1120091 

can recombine. This might be considered the beginning 
of the plateau phase. 

On the plateau, the position of the photosphere is de- 
termined by the location of the hydrogen recombination 
front, which occurs at a temperature T\ « 6000 K. As 
the ejecta cool, the front recedes inward (Figure [3]). At 
around day 120, the front reaches the base of the hy- 
drogen envelope. Recombination occurs more quickly in 
the helium layers, so the remaining internal energy is 
depleted quickly and the light curve drops off sharply, 
ending the plateau. 

The evolution of the photospheric velocity over time 
(Figure S|) a grees well with the o bservations of SNe IIP 
compiled bv lNugent et al.l (j2006). Had we ignored the 
effects of hydrogen recombination, the photospheric ve- 
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o.o I ... i ... i ... i ... i ... i ... i ... i . I 

4000 4200 4400 4600 4800 5000 5200 5400 
Wavelength (Angstroms) 

Fig. 5. — Evolution of the Fe II and Ti II lines in the spectral 
range 4000-5550 A. The recession of the photosphere is observable 
in the decreasing Doppler shift of the line minima, especially that 
of Fe II A5169 line. 

locity would have declined much more slowly, in conflict 
with the observed. Thus, although the photospheric ve- 
locity is often taken as a measure of (E /M^) 1 / 2 (i.e., 
Eq. [T]) for SNe IIP this is clearly only valid for times be- 
fore recombination sets in (here t < 20 days). At later 
times, the position of the photosphere is largely deter- 
mined by the inward progression of the recombination 
front. 

Figure [6] further illustrates how the nature of the opac- 
ity affects the light curve. If we artificially increase the 
electron scattering opacity by a factor of 2, the light curve 
becomes dimmer and broader, in agreement with the an- 
alytical scalings of Eq. [H] This indicates that radiative 
diffusion in the ionized regions is indeed significant. If 
we neglect the effects of hydrogen recombination, the re- 
sulting light curve declines in a roughly power law fash- 
ion, with a lower average luminosity and longer duration. 
This implies that the recombination wave is responsible 
for the flatness of the plateau and the steep drop off af- 
terward. 

The opacity is also affected by the helium abundance 
in the hydrogen envelope. Because helium recom- 
bines at higher temperatures than hydrogen, a larger Xh e 
effectively reduces the electron scattering opacity. The 
light curve of a model with X\ lc = 0.5 is therefore 0.4 mag 
brighter and 20 days shorter than one assuming pure hy- 
drogen (Figure [7]) ■ Helium in the core of the ejecta also 
affects the light curve, though in a slightly different way. 
For models with larger helium cores, the recombination 
front will reach the base of the hydrogen layer at an ear- 
lier time, and so the plateau will end relatively sooner. 

As expected from the analytical arguments of Sec- 
tion [2] the inclusion of radioactive 56 Ni extends the 
plateau duration, but has essentially no effect on the lu- 
minosity at times < 50 days (Figure Because 56 Ni is 
synthesized only at the ejecta center, radioactive energy 
does not have enough time to diffuse out and affect the 
plateau unless extreme masses or outward mixing of 56 Ni 
are considered. 

6.2. Broadband Light Curves and Spectra 

The model broadband light curves are shown in Fig- 
ure [3 The U and B band light curves decline sharply, 
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Fig. 6.— V-band light curve of model M15.E1.2_Z1 (with no 
56 Ni included) computed using different opacity prescriptions. The 
solid line properly includes all relevant opacities. The dash-dotted 
line shows a model in which hydrogen is not allowed to recom- 
bine, which demonstrates the importance of ionization effects on 
the plateau. The dashed line shows a model in which the electron 
scattering opacity was increased by a factor of 2, which demon- 
strates the importance of diffusion in the ionized regions. 
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Fig. 7. — Bolometric light curves of model M15_E1.2_Z1 (with a 
56 Ni mass of 0.06 Mq) computed with different helium abundances 
in the hydrogen envelope. A greater helium abundance reduces the 
electron scattering opacity and leads to a shorter, brighter plateau. 



showing virtually no plateau, while the V, R and / bands 
are flatter. This behavior, which is also seen in observa- 
tions, can be understood by examining the spectral evo- 
lution on the plateau (Figure [9|). At longer wavelengths 
(A > 5000 A) the continuum is fairly well approximated 
by a blackbody of constant temperature. The V — I and 
V—R colors are therefore fairly constant over the plateau. 
At shorter wavelengths (A < 5000 A) on the other hand, 
the spectrum is heavily affected by the blanketing from 
millions of blended iron group lines (in particular those 
of Fe II and Ti II) reflecting the metallicity of the pro- 
genitor star. This line opacity depends sensitively on 
temperature, as a slight cooling of the photosphere in- 
duces Fe III and Ti III to recombine to Fe II and Ti II. 
The corresponding non-linear increase in line blanket- 
ing, clearly visible in Figure [H causes a drop in U and B 



Type IIP Supernova Light Curves 



7 




50 100 150 200 

Days Since Explosion 



Fig. 8. — Broadband UBVRI light curves of our standard model 
M15_E1.2_Z1 (with an ejected 56 Ni mass of 0.13 M Q ). The faster 
decline in the U and B bands is due to increasing line opacity (see 
Figure [9} . 

magnitudes much greater than would be expected from 
a pure blackbody spectrum. 

Figure D_] illustrates how the inward progression of the 
supernova photosphere is detectable in the spectral se- 
ries. The Doppler shifts of Fe II and Ti II absorption 
lines in the wavelength region 4000-5000 A decrease over 
time. In most applications, the Fe II A5169 line is used 
to infer the photosphcric velocity, as it is strong enough 
to be measured relatively easily, but weak enough to not 
be saturated above the photosphere. 

7. MODEL SURVEY 

Table 2 and Figures [TO] and QT] summarize the light 
curve properties of the entire model survey. The models 
vary by more than a factor of 10 in plateau luminosity, 
and by about a factor of 2 in duration. It is immedi- 
ately clear that most of the variation in SNe IIP events 
reflects differences in explosion energy - changes in pro- 
genitor mass only account for a factor of ~ 2 in lumi- 
nosity. By directly comparing to the observed sample of 
nearby SNe IIP (see Section [9]), we infer that the explo- 
sion energy of real SNe IIP spans the range 0.6 — 4.8 B, 
with a typical mean value around 0.9 B. 

Of the analytical scaling laws discussed in Section [21 
the model luminosity dependence follows most closely 
those of Eqs.[5J and in particular L sn oc E 5 / 6 . The scaling 
of the plateau duration, however, deviates from Eqs. [5] 
following more closely t p oc E 1 / 4 . Guided by the analytic 
results, we find expressions that well fit the models 

L 50 = 1.26 x 10 42 El( G M- 1/2 Rl% X^ e ergs s" 1 , 

tp,o = 122 i^XoXSo^e 2 days, 

where i?o,500 = i?o/500i?© and t Pi o is the plateau du- 
ration when no 56 Ni is included. The dependence 
accounts for the effects of helium both in the envelope 
and the core. Figure [12] illustrates that the accuracy of 
these expressions is quite good. 

In principle, Eqs. [TTJ along with the expression for the 
scaling velocity Eq. [JJ could be used to infer the physical 
parameters (E, M, Rq) from the o bserved (L, t p , v), in the 
manner applied by Hamuy (2003). In practice, there are 
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Fig. 9.— Spectral evolution of model M15.E1.2_Z1, with the 
time since explosion marked on the figure. At longer wavelengths 
(A > 5000 A) the spectrum is well approximated by a blackbody at 
T ~ 6500 K (dotted lines). At shorter wavelengths, the iron group 
line blanketing becomes progressively stronger over time, reducing 
the flux. 
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Fig. 10. — Model V-band light curves from our survey, assum- 
ing no 56 Niejected. Top: light curves of the M = 15 Mq solar 
metallicity progenitor with different explosion energies (marked on 
the figure in B). Higher energy explosions are brighter and shorter. 
Bottom: light curves of 1.2 B explosions with different progeni- 
tor initial masses (marked on the figure in Mq). There is a non- 
monotonic behavior in the M = 25 Mq light curve (dashed line) 
due to its increased mass loss and envelope helium abundance. 
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TABLE 2 
Properties of Supernova Models 



Name 


Mi 


E 


M p 


St 


Af„; 


Mm 


^50 


tp,0 






M12_E1.2_Z1 


12 


1.21 


1 


36 


9.53 


0.16 


1.91e42 


116 


-17.25 


4915 


M12_E2.4_Z1 


12 


2.42 


1 


36 


9.53 


0.18 


3.67e42 


99 


-17.98 


6346 


M15_E1.2_Z1 


15 


1.21 


1 


48 


11.29 


0.26 


2.16e42 


124 


-17.38 


4959 


M15.E2.4_Z1 


15 


2.42 


1 


48 


11.29 


0.31 


4.35e42 


105 


-18.15 


6491 


M15.E0.6.Z1 


15 


0.66 


1 


48 


11.25 


0.24 


1.26e42 


149 


-16.79 


3966 


M15.E4.8.Z1 


15 


4.95 


1 


48 


10.78 


0.36 


7.80e42 


88 


-18.80 


8479 


M15.E0.3.Z1 


15 


0.33 


1 


48 


11.27 


0.22 


5.93e41 


177 


-15.96 


3125 


M20.E1.2_Z1 


20 


1.22 


1 


54 


14.36 


0.34 


2.61e42 


111 


-17.57 


4947 


M20.E2.4_Z1 


20 


2.42 


1 


54 


14.37 


0.40 


4.85e42 


119 


-18.26 


6459 


M20.E0.6.Z1 


20 


0.68 


1 


54 


14.36 


0.32 


1.40e42 


167 


-16.89 


3979 


M20.E4.8.Z1 


20 


4.99 


1 


54 


14.37 


0.48 


8.57e42 


99 


-18.91 


8337 


M25.E1.2_Z1 


25 


1.22 


1 


59 


14.22 


0.37 


3.94e42 


131 


-18.00 


5033 


M25_E2.4_Z1 


25 


2.43 


1 


59 


14.22 


0.43 


6.66e42 


107 


-18.59 


6483 


M25.E0.6.Z1 


25 


0.66 


1 


59 


14.11 


0.34 


1.96e42 


154 


-17.23 


4281 


M25.E4.8.Z1 


25 


5.00 


1 


59 


12.97 


0.56 


1.10e43 


86 


-19.17 


7948 


M15_E1.2_Z0.1 


15 


1.26 


1 


16 


13.27 


0.12 


1.67e42 


130 


-17.04 


4716 


M15.E2.4.Z0.1 


15 


2.48 


1 


16 


13.24 


0.16 


3.08e42 


107 


-17.71 


6098 


M15.E0.6.Z0.1 


15 


0.65 


1 


46 


13.28 


0.10 


8.59e41 


156 


-16.32 


3671 


M15.E4.8.Z0.1 


15 


4.90 


1 


16 


13.18 


0.20 


5.31e42 


88 


-18.30 


7670 



several complicating factors. The envelope helium abun- 
dance and the size of the helium core, for instance, are 
significant factors, but unfortunately there are no clean 
observables to constrain them. In addition, the photo- 
spheric velocity on the plateau is largely determined by 
recombination, and thus not necessarily a good measure 
of the ejecta velocity v sn oc (E/M^) 1 / 2 (see Figure|4|). A 
measurement of the velocity at epochs prior to recombi- 
nation (t < 20 days) is therefore preferred. 

An alternative approach would use the fact that in 
the progenitor models, i?o, M e j, and Xn c are correlated, 
so that some of the degeneracies may be removed. For 
future photometric surveys, useful relations would allow 
for a determination of E and Mj given only L50 and i Pi o. 
We find for solar metallicity models 

L 50 = 1.49 x 10 42 E°f 2 M^ ergs 

t P ,o — 128 E~°' 26 M^i Q days, 

where Mi nj io = Mi/10M Q . These relations (which fit the 
models to within 10%) can be applied to infer the gross 
properties of SNe IIP without need for follow-up spec- 
troscopy. However, one should bear in mind that they 
rely on the predictions of stellar evolution and explo- 
sion calculations, and thus are subject to uncertainties 
in, e.g., mass loss and fallback. 

Before applying either Eqs.[TT]or Eqs. [12] it is critical 
to account for the fact that 56 Ni in the ejecta tends to 
extend the plateau. Figure [13] shows that our derived 
analytical scaling (Eq. [TCT| fits reasonably well, with the 
refined numerical values 

/ \ 1/6 

t P = t pfi x II OMMwE-^R^M^ 2 J . (13) 

We find that the luminosity on the tail of the light curve 
is nearly identical to the instantaneous energy deposition 
from 56 Co decay. The ejected mass of 56 Ni can then be 
inferred in the typical way, by measuring the luminos- 
ity at a point on the tail. Unfortunately, M c j and Rq 
also appear in this expression; however, their approxi- 
mate values for a given initial mass could be taken from 
Table 1. 
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Fig. 11. — Light curve properties of all of the models in our 
survey (assuming no 56 Ni ejected). Left: bolometric luminosity on 
the plateau (measured 50 days after explosion) as a function of the 
explosion energy. Circles denote solar metallicity models, squares 
0.1 solar models, and the size of the symbol is proportional to 
the progenitor star initial mass. The solid line shows the L oc _E 5 / 6 
scaling. Right: same as the left panel, but for the plateau duration. 
The solid line shows the t p oc i? 1 / 6 scaling, the dashed line the 
t p oc E 1 / 4 scaling. 




Fig. 12. — Accuracy of the analytic scalings Eqs,l8lfor the plateau 
luminosity (left) and duration (right). The x-axis values are those 
determined by the numerical simulation, while the j/-axis values 
are those derived from the analytic equation. Circles denote solar 
metallicity models, squares 0.1 solar models, and the size of the 
symbol is proportional to the progenitor star initial mass. 

8. BOLOMETRIC CORRECTIONS AND DUST 

From the models, one can derive formulae useful for 
making bolometric and dust corrections to observations. 
Figure [14] plots the difference in bolometric and F-band 
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Fig. 13. — Effect of radioactive Ni on extending the plateau 
duration. The j/-axis is the plateau duration t p divided by the t p 
for the zero 56 Ni case. The models included here have Mi = 15 Mq, 
explosion energies of 1.2 and 2.4 B, and 56 Ni mass varied from 
to 0.35 Mq. The solid line shows the analytic scaling Eg. 1131 



Fig. 14. — Bolomctric corrections (i.e., bolometric magnitude 
minus V-band magnitude) at day 50 for all models. Circles are 
solar metallicity models and squares Z = 0.1 solar model. The 
dashed line shows a linear fit to the solar metallicity models, Eg. 1141 



magnitude at day 50 for all models. The typical bolo- 
metric corrections are around 0.2 mag, but increase for 
brighter events by as much as 0.07 mag. For solar metal- 
licity models, we fit the relation 



BC 50 = 0.24 - 0.02 5 x (M v ,so + 18), 



(14) 



where V50 is the V-band magnitude at day 50. The 
Z = 0.1 solar metallicity models, due to the lesser line 
blocking, have bolometric corrections about 0.07 mag 
lower. 

Previous studies have typically estimated dust extinc- 
tion by measuring the V — I excess over an assumed in- 
trinsic color. We find the models have a roughly constant 
color on the plateau of V — I « 0.5; however, there is a 
slight trend for brighter models to be bluer (Figure [T5|) . 
A fit to the solar metallicity models at day 50 gives 

(V - 7) 5 o = 0.52 + 0.03 x (V 5Q + 17.5). (15) 

The values are simi lar to the value (V — I)o = 0.53 that 
iNueent et all (j2006fl inferred by examining the ridge line 
of observed events. 

As mentioned in Section [5j the model /-band mag- 
nitudes are sensitive to NLTE effects, especially in the 
Ca II IR triplet line. Eg. [TBI was therefore determined us- 
ing day 50 spectrum calculations which treated calcium 
in NLTE, though under the stationary approximation. 
The model predictions are thus subject to uncertainties 
in the assumed calcium a bundance and perhaps to time 
dependent NLTE effects (jDessart fc Hillier|[200l . 

9. STANDARD CANDLE RELATIONSHIP 

Using a sample of nearby observed SN IIP, 
lHamuv fc Pinto! (|2002f ) found that the plateau luminos- 
ity (measured at day 50 after the explosion) correlated 
rather tightly with the photospheric velocity, as mea- 
sured from the Doppler shift of spectral absorption lines. 
This empirical standard candle (SC) relation provides a 
simple means for calibrating SN IIP luminosity for dis- 
tance measures. 

Figure [TBI shows the Hamuy SC relation for our model 
survey set, here in terms of the V-band magnitude (My) 
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Fig. 15. — The V — I color at day 50 for all models. Circles 
denote solar metallicity models and squares 0.1 solar metallicity 
models. Th e so lid line shows a linear fit to the solar metallicity 
models, Eq. [15] 

and the photospheric velocity. The model relation is as 
tight or tighter than the observed, and with a similar 
slope. The rms dispersion is only a = 0.27, which trans- 
lates to ~ 13% errors in distance measures. To first or- 
der, the velocity and luminosity of SNe IIP are both set 
by the explosion energy. The dispersion in the relation is 
due to variations in the progenitor mass and metallicity 
for a given explosion energy. 

The physical interpretation of the model SC rela- 
tion is straightforward, being essentially a recasting of 
the Baade-Wesselink or expanding photosphere methods 
that have been in use for many years. The luminosity is 
written using Stefan's law and the radius of the super- 
nova photosphere R = v p ht 



1 ph> 



(16) 



where C is a "dilution" factor which accounts for devia- 
tion of the spectrum from blackbody. In Type II atmo- 
spheres, both t he effects of scattering and line blank eting 
contribute to £ (|Wagonerlll98H[Eastman et ail 19961 ). To 
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Fig. 16. — Model standard candle relationship between the V-band plateau luminosity (measured at 50 days) and the photospheric 
velocity. Left: the models are color coded by explosion energy. Circles denote solar metallicity models, squares 0.1 solar models, and the 
size of the symbol i s proportiona l to the progenitor star mass. Right: Comparison of the model relation (circles) with the observation 
sample compiled by ( Hamuy 2003) (open squares). The large filled square is SN 1999em using the Cepheid distance of Leonard ct al. (2003) 
and an extinction of A v = 0.31. 



determine L using the expanding photosphere method, 
the observer measures i> p h and the time since explosion t, 
and estimates the photospheric temperature T p h from the 
color of the spectrum. The dilution factor must be cal- 
culated using detailed numerical models (the main com- 
plexity of the approach) . NLTE spectral modeling finds 
that £ varies between 0.5 and 2.0, and is chiefly a func- 
tion of luminosity, being rather insensitiv e to other ejecta 
parameters such as the density structure (|Eastman et al.l 
1996; IDessart fc Hiilienl2005al) . 

The standard candle relation is simply an expression 
of Eq. [16] under certain restricted conditions. The time 
since explosion t is, by construction, fixed at 50 days. 
The temperature T pn for SNe IIP on the plateau is nearly 
a constant, constrained to be near the recombination 
temperature Ti k, 6000 K. The dilution factor £ may 
vary from event to event, but if £ is primarily a function 
of luminosity this dependence can be absorbed into the 
exponent. This implies L — Cv"~^ e , where the constant 
C and the non-blackbody effects e can be calibrated us- 
ing a sample of nearby objects, or a set of theoretical 
models. 

The SC relation need not be applied only at day 50, 
and we find that similar relations apply all along the 
plateau. However the time since explosion must be 
known as the normalization depends on time (Eq. I16j) . 
We find that an uncertainty in explosion time of 10 days 
leads to an error in inferred brightness of 0.2 — 0.3 mag. 
It is unwise to apply the SC relation at times much ear- 
lier than 30 days, as the ejecta temperatures are likely 
too high for recombination to have set in, and there is 
no assurance that T p h ss Ti. 

One nice feature of the models is that they offer an ab- 
solute normalization of the SC relation without needing 



to assume a value of the Hubble constant. By fitting the 
relation evaluated at different times since explosion, we 
find 



M v {t) = -17.4- 6.9 log 10 (v ph (t)/5000.0) 
+3.11og 10 (t/50 days) 



(17) 



The models do predict a deviation from the simple 
L oc Up h relation of Eq. [16J showing instead L oc v^ 5 in 
general accordance with that f ound in the observational 
sample ([Hamuv fc P into 2002). This effect is primarily 
due to the deviation of the spectrum from a blackbody. 

The model relation of Figure [16j has a similar n orma l- 
ization to the observations, taken from iHamuvl (2003). 
This implies that our model SC relation is in rough 
agreement with the distances to SNe IIP obtained in 
other ways. Particularly comforting is the agreement 
with SN 1999em, which has a measured Cepheid dis- 
tance to its host gal axy NGC 1637 of 11.7 ± 1.0 Mpc 
([Leonard e t al. 2003|). We find a very similar distance 
of 11.6 ± 1.2 Mpc from Eq. [17] when taking the ob- 
served va lues m v = 13.98. y n h = 3757 km s -1 , and 
(following iBaron et all feOOO: Ha muv et aTll2001h . an ex- 
tinction of A v = 0.31. This distance is also consistent 
with independent estimates using the exp anding pho- 
tosphere method (IDessart fc Hillieil l2005al) and SEAM 
(jBaron et al.|[200l . 

One drawback of the standard candle method, from 
the observational point of view, is that a high quality 
spectrum is needed to measure the photospheric velocity 
- a difficult prospect for high redshift events. As future 
surveys will observe light curves for a enormous number 
of SNe IIP with limited spectroscopic follow-up, meth- 
ods of purely photometric calibration, however coarse, 
may be of interest. As the explosion energy is the pri- 
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Fig. 17. — Relationship between the plateau duration (assum- 
ing zero 56 Ni) and the luminosity at day 50. Circles denote solar 
metallicity models, squares 0.1 solar models, and the size of the 
symbol is proportional to the progenitor star mass. 

mary variable determining both the plateau luminosity 
and duration, we explored the relationship between these 
two observables. A relationship exists (Figure fT7|) and is 
fit by 

M v ,50 = -18.4 - 0.03[i p ,o - 100] (18) 

Applying this relation reduces the dispersion from 1 mag 
down to 0.4 mag. In practice, the measured plateau du- 
ration t p must be corrected for the effect of the ejected 
56 Ni mass on its duration in order to determine i p .o- The 
residual scatter in the relation is clearly due to variation 
in progenitor initial mass or metallicity for a given ex- 
plosion energy. Presumably, the scatter could be reduced 
further by using additional light curve relation, such as 
the color evolution. 

10. DISCUSSION AND CONCLUSIONS 

We explored the light curves and spectra of SNe II 
models with various progenitor masses, metallicities, and 
explosion energies. We found that explosions with ener- 
gies 0.3 — 4.8 B of stars with initial masses in the range 
12 — 25 Mq can explain the observed range of luminosi- 
ties, velocities, and light curve durations of most SNe IIP. 
For existing and future observational surveys, the model 
results should be useful for inferring the progenitor star 
properties, explosion energies, distances, and dust ex- 
tinction of observed events. 

This study, as have previous studies, quantified how 
the basic supernova parameters (M e j, Rq, and E) affect 
the light curves. We also highlighted the important role 
of two additional parameters: the radioactive 56 Ni mass 
and the envelope helium abundance. The presence of 
56 Ni extends the plateau duration, but typically does not 
affect the luminosity on the plateau at times t < 50 days. 
The neglect of the e ffect of 56 Ni may be the main reason 
why lHamuvl {2003), in his analysis of 16 SNe IIP, in- 
ferred implausibly large ejecta masses (up to 50 M ). In 
that study, the longer plateau duration would have to be 



accounted for by an increased diffusion time, and hence 
larger ejecta mass. Here we presented analytical formu- 
lae which may be useful in accounting for the effects of 
56 Ni on the plateau. 

The models confirm the standard candle method of 
calibrating SNe IIP and illuminate its physical origin. 
The method is a promising way to determine distances 
to SNe IIP, with a clear physical explanation in terms of 
the ionization physics of hydrogen. On the other hand, 
the models raise some concerns about systematic errors. 
Progenitors with different masses or metallicities lie on 
differently normalized relations in Figure 1161 If the pro- 
genitor population at high redshift has different demo- 
graphics than that at low redshift (as might be expected) 
a systematic bias may be introduced into distance mea- 
surements. The effect of going from Z = 1 to Z = 0.1 
solar metallicity in the models is at the 0.1 mag level. 
It may be possible to reduce these errors by using color 
information from the light curve. 

We find a correlation between plateau luminosity and 
plateau duration which could be useful in roughly cali- 
brating SNe IIP luminosities using only photometric data 
(to about 20% in distance) . This correlation reflects the 
fact that in the models one parameter, the explosion en- 
ergy, primarily controls both the light curve brightness 
and duration, while the progenitor star properties play a 
secondary role. The validity of such a relation needs to 
be empirically checked, as the scatter will be smaller or 
larger depending on whether the bulk of SNe IIP arise 
from a narrower or wider range of progenitor masses and 
radii than that considered here. In practice, the rela- 
tion also needs to take into account the effect of 56 Ni on 
extending the plateau duration. 

Correction for dust extinction remains a difficult issue 
for determining the distances to SNe IIP. The models 
provide some theoretical guidance as to the intrinsic color 
evolution of SNe-IIP light curves, however their accu- 
racy may be limited by the assumptions in the radiative 
transfer, and are sensitive to variations in the envelope 
metallicity. On the other hand, one could try to invert 
the problem. Assuming the cosmological parameter are 
accurately constrained by other means, one could use the 
standard candle method to solve for the dust extinction 
of SNe IIP, thus providing an estimate of the variation of 
dust properties with galactic environment and redshift. 
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